Some remarks on Kurepa's left factorial 

Bernd C. Kellner 



Abstract 

We establish a connection between the subfactorial function S(n) and the left fac- 
torial function of Kurepa K(n). Some elementary properties and congruences of both 
functions are described. Finally, we give a calculated distribution of primes below 
10000 of K(n). 
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1 Introduction 

The subfactorial function is defined by 
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which gives the number of permutations of n elements without any fixpoints, also 
called derangements of n elements, see [6, p. 195]. This was already proven by P. R. 
de Montmort [2] in 1713. L. Euler [3] independently gave a proof in 1753, see also [4]. 
This function has the properties (e ~ 2.71828 is Euler's number) 
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Kurepa's left factorial function is defined by 



n-l 



fc=0 



In 1971 D. Kurepa [8] introduced the left factorial function which is denoted by \n = 
K{n). Sometimes the subfactorial function is also denoted by !n, so we do not use this 
notation to avoid confusion. For more details of the following conjecture see a overview 
of A. Ivic and Z. Mijajlovic [7]. 



Conjecture 1.1 (Kurepa's left factorial hypothesis) 

The following equivalent statements hold 



(K(n),n\) = 2, n>2, 

K{n) ^ (mod n) , n>2, (KH) 
K (p) ^ (mod p) , p odd prime. 



Recently, D. Barsky and B. Benzaghou [1] have given a proof of this hypothesis. Since 
K(n) is also related to Bell numbers B n via 

K(p) = Bp-\ — 1 (mod p) 

for any prime p, they actually proved that B v -\ ^ 1 (mod p) is always valid for any 
odd prime p. 

Gessel [5, Sect. 7/10] gives some recursive identities of S(n), B n , and others with umbral 
calculus. Define symbolically S n = S(n) and B n = B n with S° = B° = 1, then one 
may write 

B n+1 = (B + l) n and n\ = {S + l) n , n>0. (1.4) 
Interestingly, both sequences have the same property as follows. 
Lemma 1.2 Let p be a prime. Then 

]T(-l) fc £ fc ^]T(-l) fc S(fc)^0 (modp) 

k=0 k=0 

with 

Bp = 2 (mod p) and S(p) = —1 (mod p) . 



Proof. By (1.1) and Wilson's theorem, we have S(p) = —1 = (p — 1)1 (mod p). 
Hence, we can rewrite (1.4) by B p = (B + l) p ~ l and S p = (S + l) p ~ l (mod p). Since 
( p ^ 1 ) = (— l) k (mod p) for < k < p, this provides the proposed congruence. Now, we 
use a congruence of Touchard for Bell numbers, see [5, Sect. 10, Theorem 10.1]. Then 

B n+P - B n+ i - B n = (mod p) , n>0. 

With n = and Bo = B\ = 1, we obtain B p = 2 (mod p). □ 

First values of K(n), S(n), and B n are given in the following table. 
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2 Congruences between K(ri) and S(ri) 

Lemma 2.1 Let n be a positive integer, then 

K(n) = (-l^Sin - 1) (mod n) . 

Proof. Case n = 1 is trivial. Let n > 2. Then we have 

n— 1 / 1N n— l 



(_i r -l 5(n _l )= ^(_i r -l- 



n — l 
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l-fc)! = £(-!) 



re — l 



fc=o v ' fc=o 

by turning the summation. Since it is valid for < < re 

(-l) fc f n ~ k ^jk\ = {-lf{n _ l) ...(„- fc) = fc! (mod n) 

this provides, term by term, the congruence claimed above. 



k\ 



□ 



By Lemma 2.1 and (1.3), we easily obtain a generalization, however, which is only 
noted for primes elsewhere. 



Corollary 2.2 Let n be a positive integer, then 

(n-l)! 



n-l 



+ 5 n -i (mod n) 



tf(re) — . (-1) 

Hence, (KH) is equivalent to 

(re - l)! 

^ — (mod re) , re > 2 , 

e 

while by recursive property (l.l) 

re! 

— = — S n -\ (mod re) , re > l 

e 

is always valid. 

Corollary 2.3 Let n be a positive integer, then (KH) is equivalent to 

(mod re) re =1,2. 
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Lemma 2.4 Le£ p be a prime, then 



K{p)-K(p-l) = -^-^ (modp), Z = l,...,p. 
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PROOF. Let I e {1, . . . ,p}. We then have 

/ 



p-l 
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K(p)-K(p-l)= fc!=&-fcNE([4! = -f^r 



k=l 
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since 



(p-k)\ 



(fc-l)! 



(mod p) 



follows by Wilson's theorem. 



(2.1) 



□ 



Corollary 2.5 Let p be an odd prime, then (KB.) implies for < / < p 



K(p-l-l)£ 



S(l) 



(mod p) 



respectively 



l\K(p-l-l) # 



+ 5[ (mod p) . 



Since (KH) is true, we obtain, as an example, the following congruences 
K(p)^0, K(p-l)^l, tf(p-2)#0, K(p-3)^±, K(p-4)^± (mod p) . 



3 Properties of K(n 



To describe some interesting properties of K(n), we introduce the following definition 
which we name after Kurepa. 

Definition 3.1 Let p be an odd prime. The pair (p,n) is called a Kurepa pair if 
p r | K(n) with some integer r > 1. The max. integer r is called the order of (p,n). 
The index of p is defined by 

*if(p) = #{ n : {P-, n ) is a Kurepa pair} . 
If ik(p) > 0, then p is called a Kurepa prime. 

We have, e.g., the Kurepa pairs (19, 7), (19, 12), and (19, 16). If (KH) would fail at an 
odd prime p, then this would imply ixip) = oo. This is an easy consequence of 

p | K(p) , p | (p + m)\ for m > . 

The case p = 2 is handled separately. One easily sees that 2 | K(n) for n > 2 and 
If (re) = 2 (mod 4) for n > 4. First values of ix{p) are given in the following table. 
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Theorem 3.2 Let (p, n) be a Kurepa pair. Then p > n > 3 is valid with 

K{p) = {-l) n n\S{p-l-n) (mod p) 

which implies p \ S(p — 1 — n). Furthermore one has ixip) < L^r^J- Consequently, 
there exist infinitely many Kurepa primes. 

Proof. For now, let p be an odd prime. Let (p, n) be a Kurepa pair. Since p \ K(p+m) 
for m > by validity of (KH) and first values of K(-) are 0,1,2,4, this yields p > n > 3. 
We use Lemma 2.4 with n = p — I, then we have 

£ K(p) = K(p) - K{n) = - S(P ~ 1 ~ (modp) 

(p — 1 — raj ! 

which provides the result by means of (2.1) and also p f 5(p — 1 — n). Now, we have 
to count possible Kurepa pairs. Corollary 2.5 shows that K{p — 2) ^ (mod p). If 
p | K(n) then p { ET(n + I) for / = 1, 2, 3. This is seen by n! ^ (mod p) and 

n\ + (n + 1)! = (n + 2) n! ^ , n! + (n + 1)! + (n + 2)! = (n + 2) 2 n! ^ (mod p) , 

since n ^ p — 2. On the other side, we have 4 < n < p — 1. Then a simple counting 
argument provides ik(p) < L^i^"J- Finally, K(n) — > oo for n — > oo and p | K(n) =^> 
p > n for odd primes imply infinitely many Kurepa primes. □ 

Now, the remarkable fact of K(n) is the finiteness of Kurepa pairs for all odd primes. 
In p-adic analysis, the series 

oo 

K{oo) = J2k\ 

k=0 

is an example of a convergent series resp. K{n) is a convergent sequence which lies in 
Z p , the ring of p-adic integers. Then (KH) is equivalent to K{oo) is a unit in Z p for 
all odd primes p. The behavior (mod p r ) is illustrated by the following theorem. Note 
that l r is related to the so-called Smarandache function for factorials. 

Theorem 3.3 Let p, r be positive integers with p prime. Then the sequence 

K{n) (mod p r ) , n > 

is constant for n > l r p with r > l r and 

L = mm < / : / H > r 

I { p- 1 

where o~ p {l) gives the sum of digits of I in base p. 

Proof. We have to determine a minimal I with the property ordj, (Zp)! > r. Counting 
factors which are divisible by p, we obtain 

ord p (Zp)! = Z + ord p Z! = Z + - - ^ 

by means of the p-adic valuation of factorials, see [9, Section 3.1, p. 241]. □ 
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At the end, we give some results of calculated Kurepa pairs. There are N = -/r(lOOOO) — 
1 = 1228 odd primes below 10000. Let N r be the number of odd primes with index 
ik(p) = r in this range. The following table shows the distribution of the index ix- 



r 





1 


2 


3 


4 


5 


N r 


459 


472 


213 


58 


23 


3 


N r /N 


0.37378 


0.38436 


0.17345 


0.04723 


0.01873 


0.00244 



The calculated Kurepa pairs with index ix(p) = 5 are as follows. 



(2203,277) 


(2203,788) 


(2203,837) 


(2203,1246) 


(2203,1927) 


(5227,850) 


(5227,1752) 


(5227,3451) 


(5227,4363) 


(5227,4716) 


(6689,1716) 


(6689,2404) 


(6689,3641) 


(6689,3969) 


(6689,6601) 



All primes below 10000 appear with a simple power in K(n), except K(3) = 4. On 
the other side, the occurrence of higher powers p r in K(n) seems to be very rare. 
M. Zivkovic [10] gives the first example 54503 2 | if(26541). There are two Kurepa 
pairs (54503,26541) and (54503,49783), but only the first of them has order two. 

One may ask whether the distribution of Kurepa pairs resp. the index ix can be 
asymptotically determined and even proven. Are there infinitely many non-Kurepa 
primes p with ix(p) = 0? It seems that this subject of K(n) and its distribution of 
primes will be much simpler to attack as, for example, the more complicated but in a 
sense similar case of the distribution of irregular primes of Bernoulli numbers. 
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